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Introducción 





Con el objetivo de poner de manifiesto ciertos rasgos 
de la imagen que pueden quedar enmascarados con 
la abundancia de datos, o bien con el fin de destacar 
algún otro tipo de información poco evidente o in- 
cluso aparentemente oculta, se pueden realizar ope- 
raciones algebraicas entre bandas de una misma ima- 
gen o entre imágenes diferentes. Mediante el álgebra 
de bandas o de imágenes se sintetiza en un solo valor 
información procedente de varias o de todas las ban- 
das de la imagen. 


Las transformaciones globales, es decir, aquellas en 
las que toda la imagen resultante depende de toda la 


imagen o imágenes originales, se basan en el álgebra. 
En general, el álgebra de imágenes proporciona nue- 
vas imágenes con bandas obtenidas a partir combina- 
ciones lineales de las bandas originales. 


Existen multitud de operaciones que pueden reali- 
zarse con imágenes o con bandas de ellas, orientadas 
cada una a un objetivo concreto. En esta lección se 
estudiarán el análisis de componentes principales, la 
transformación Tasseled Cap, los cocientes de ban- 
das y, como caso particular de ellos, los índices de 
vegetación. 


Análisis de componentes principales 





El análisis de componentes principales (ACP) o trans- 
formación a componentes principales es una técnica 
estadística multivariante muy utilizada para definir la 
mejor combinación de nuevas bandas capaz de con- 
servar la mayor parte de la información de la imagen 
completa. 


El ACP busca los ejes de mayor variabilidad e infor- 
mación de la imagen y es de gran utilidad para los 
estudios multitemporales, o para seleccionar tres 
nuevas bandas para una composición en color que 
concentren más información que cualesquiera de las 
originales. 


La distribución de niveles digitales en una banda de 
una imagen multiespectral puede considerarse una 


variable aleatoria y, como tal, puede someterse a cual- 
quiera de las técnicas estadísticas habituales 


La función de densidad univariante de la distribución 
de los ND es en realidad el histograma de cada banda. 
Para facilitar el trabajo con los datos, dicha distribu- 
ción pue asimilarse a una distribución normal. De 
esta forma, la distribución de los ND cada banda 
queda caracterizada completamente por sólo dos pa- 
rámetros: la media y la varianza de sus datos. 


A su vez, considerada la imagen multiespectral com- 
pleta, los ND multiespectrales pueden considerarse 
como una variable aleatoria vectorial multivariante y 
cada banda constituye una componente del vector. 
La función de densidad multivariante es ahora el his- 
tograma multiespectral de la imagen y el equivalente 
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multiespectral del valor promedio es el vector de me- 
dias. 


La relación entre bandas está caracterizada por las 
covarianzas o los correspondientes coeficientes de 
correlación. 


En una imagen multiespectral, el vector 
laz) (1) 


Se denomina vector de características y representa 
de forma vectorial los ND de las n bandas de la ima- 
gen. 


Dentro de una determinada banda K, el valor medio 
es: 


Zi 
y, = YE (2) 


1j 


en donde N es el número de celdas de la banda, igual 
al producto del número de columnas por el de líneas. 


La varianza de esa banda es: 


(sy) m1) 

O A 

o, = 2, Ñ (3) 

La covarianza entre dos bandas K y L es: 

dui ed sl IED —m, | (4) 
ij 


. 2 
Y conviene recordar que O, ,=0,, Y que 0, =0;, 
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esto es, la varianza de la propia banda. 


Finalmente, el coeficiente de correlación entre esas 
dos bandas es: 


Se Se Se 


ko 7 Sa 
ida y/o ¿a MO 


La varianza de una banda es una medida de la dis- 
persión de sus ND en torno al valor medio y repre- 
senta el contraste de la imagen en esa banda. Visto 
desde otro punto de vista, es una buena medida de la 
cantidad de información que proporciona esa banda. 
Por su parte, la covarianza y el coeficiente de corre- 
lación entre dos bandas indican el grado de relación 
que hay entre ambas, este último de forma normali- 
zada entre -1 y 1. Las bandas correladas muestran 
distribuciones de sus hiperelipsoides de puntos en el 
espacio de características alargadas. En las incorrela- 
das, sin embargo, las nubes presentan simetría con 
respecto a ejes paralelos a los de coordenadas, aun- 
que no necesariamente tienen que tener forma esfé- 
rica O hiperesférica. 


(5) 





p 


Considerando las n bandas de la imagen, se pueden 
definir las matrices de covarianza (6) y de correlación 
(7), ambas de dimensión n x n y simétricas: 


2 
o 1 61, O O, 
2 
o ¡o o 
Y = 21 2 2n (6) 
2 
Ss O 2 ss O ñ 


banda j 


Correlación positiva 





Bandas Incorreladas 





banda li 


Correlación negativa 





Correlación entre bandas. 


1 Pr Pin 
1 
R= Pa Pon (7) 
Pr Pr E l 


La matriz de covarianzas de un conjunto de bandas 
linealmente independientes es una matriz diagonal, 
porque todas las covarianzas son nulas. 


Por lo tanto, si se quiere obtener un conjunto de ban- 
das linealmente independientes a partir de otro con- 
junto cuyas bandas no lo son, hay que buscar una 
transformación que consiga diagonalizar la matriz de 
covarianzas de este último. 


Diagonalización de la matriz de covarianzas 


Dada una matriz de transformación A, se denominan 
vectores propios o autovectores de A a los que son 
colineales con sus transformados: 


AX=2X (8) 
También se puede escribir: 
AX =MX (9) 


Siendo I la matriz identidad: 


1300 az 0 
OL su 0 
E A (10) 
0.0 1 


La expresión (9) se puede poner de la forma: 
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(A-11)X =0 (11) 
Y para que esta ecuación se verifique, se requiere re- 
solver la ecuación característica que proporciona los 
valores propios de la matriz A igualando a cero el 
determinante: 


9(A)= [A-2I|=0 (12) 


Particularizando para la matriz de covarianzas, la 
ecuación característica sería: 


$0.) =|E-21|=0 (13) 


Que, desarrollada, quedaría: 


2 
1 Ór ++. O, A 0 
2 
0, SU 2. O 0 
¡E-2]|= 21 2 2n Ñ - 
2 
Os De e O, A? 
2 
o A S1 S1, 
2 
al a SA 52, =0 
2h 
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(14) 


Resolviendo la ecuación característica se obtienen n 
soluciones 21, A2,"-", An, que son los valores propios o 
autovalores de la citada matriz, los cuales habitual- 
mente se ordenan de forma que 41 > A2 >---> An. 


Los vectores propios se obtienen resolviendo para 
cada valor propio la expresión: 
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EX, =1,X, (15) 
es decir: 
2 
o 1 51) 1, Xi 2, X; 
2 
Sn Le Sn | _ A (16) 
a Xx A 
Si S.2 o n m 1 Xi 


Y la transformación que diagonaliza la matriz de co- 
varianzas es: 


PUEP=A (17) 
donde P es una matriz de paso constituida por los 


vectores propios puestos por columnas en el mismo 
orden que los valores propios asociados a cada uno 


de ellos: 
| 

P- 2 A m2 (18) 
X Xp Xan 


y A es la matriz diagonal conformada con los valores 
propios ordenados: 


0 0 
0.) 0 

A= l (19) 
0.0 A 
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La demostración de (17) es trivial, pero queda fuera 
del alcance de este texto. 


Es interesante y, en el caso del análisis digital de imá- 
genes, muy importante considerar que la transforma- 
ción a componentes principales es una transforma- 
ción ortogonal, es decir, la traza de * es igual a la de 
A: 


NN (20) 
i=l i=l 


lo que quiere decir que la varianza, es decir, la infor- 
mación, se redistribuye sin pérdidas entre los compo- 
nentes principales. Además, a la primera componente 
le corresponde la mayor parte de la varianza de la 
imagen. De la varianza restante, su mayor parte 
queda explicada por la siguiente componente, y así 
sucesivamente. Dicho de otro modo, las nuevas ban- 
das, llamadas ahora componentes principales contie- 
nen toda la información original, pero sin redundan- 
cia entre ellas, y están ordenadas de mayor a menor 
grado de información. 


Utilidad en teledetección 


El ACP es utilizado en el análisis digital de imágnes 
para mejorar el análisis visual al introducir práctica- 
mente toda la información de la imagen en la compo- 
sición en color, que solamente puede hacer uso de 
tres bandas, si bien hay que tener en cuenta que el 
peso de información que se le da a cada canal es di- 
ferente, puesto los componentes principales retienen 
cantidades diferentes de varianza. 


En teledetección multitemporal es muy útil para es- 
tudiar de series temporales. El procedimiento que 
suele seguirse es el siguiente. Se extrae el primer 
componente principal de las bandas de la imagen de 











Cambio de base en la transformación ACP. Dos 
bandas correladas dejan de estarlo si se las so- 
mete a un giro en el espacio de características 
de modo que los nuevos ejes, componentes prin- 
cipales, se orienten según los vectores propios 
de la distribución original. 


Infrarrojo próximo 





línea de suelos 


suelo húmedo 








Rojo 


Evolución espectral de un 
suelo de cultivo 


cada una de las fechas. Posteriormente se realiza una 
clasificación con los CP. Quedarán al descubierto las 
áreas de la imagen que muestren diferencias entre 
unos componentes y otros, o lo que es lo mismo, evo- 
lución entre fechas diferentes. 


Para la detección de cambios entre dos fechas, se 
transforman a CP el conjunto de todas las bandas de 
ambas fechas. Las primeras CP explicarán la mayor 


Transformación tasseled cap 
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parte de la varianza de la imagen, y por tanto inclui- 
rán todos los aspectos comunes a las dos fechas. Las 
componentes secundarias contendrán información 
sobre aquellos fenómenos que no puedan explicarse 
en las dos fechas simultáneamente, por lo que indica- 
rán la evolución de los mismos. Las últimas compo- 
nentes principales retienen la información espuria de 
la imagen, esto es, el ruido. 





Se denomina diagrama de dispersión a la representa- 
ción de la nube de puntos espectral, esto es, de nive- 
les digitales o de las magnitudes físicas de reflectan- 
cia o radiancia, en un espacio espectral o de caracte- 
rísticas de dos o más dimensiones. Lógicamente, so- 
bre un papel solamente pueden representarse diagra- 
mas de dispersión bidimensionales, pero un sistema 
informático no tiene el menor problema para proce- 
sar diagramas n-dimensionales. 


Evolución espectral de los terrenos 
cultivados 


Se puede estudiar la posición espectral de las celdas 
de terreno cultivado en un diagrama de dispersión 
desde el estado de nascencia hasta la madurez y el 
agostamiento del cultivo. A lo largo de su período 
vegetativo el rastro que va dejando la nube de puntos 
en su evolución espectral en un espacio tridimensio- 
nal adquiere la forma de un gorro borlado (tasseled 
cap). Si solo se consideran las bandas del infrarrojo 
próximo y del rojo espectral, la base del gorro está 
definida por la línea de suelos, a lo largo de la cual 
se distribuyen los distintos tipos de suelos, de 
acuerdo a su estado (humedad e iluminación, funda- 
mentalmente). 


A medida que el cultivo se desarrolla y llega a su 
pleno vigor, los puntos tienden a aproximarse al eje 
del infrarrojo, alejándose de la línea de suelos. 
Cuando se inicia el período de senescencia, la nube 
de puntos va convergiendo hacia un punto único, que 
alcanza cuando el cultivo llega al estado de marchitez 
permanente, y cuya posición tiende a regresar de 
nuevo a la línea de suelos. 


La transformación tasseled cap (TCT, Tasseled Cap 
Transformation) es una herramienta que permite 
comprimir la información espectral de una imagen en 
pocas bandas asociadas a ciertas características de la 
escena, sin perder información significativa. Consiste 
en un giro en el espacio espectral mediante el cual se 
pone de manifiesto la evolución espectral de los cul- 
tivos anuales o bienales, desde su nascencia hasta la 
madurez, el agostamiento y la marchitez final. 


Tras el giro, los nuevos ejes en un espacio espectral 
tridimensional representan las variables brillo, ver- 
dor y humedad. Los ejes brillo y verdor definen el 
plano de la vegetación y los ejes brillo y humedad, el 
de los suelos. 
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Esta transformación es empírica y se aplica con dife- 
rentes matrices a cada imagen de satélite, depen- 
diendo del número de bandas y de la longitud de 


onda de cada una de ellas. En los ejemplos siguientes 
se muestran las transformaciones sobre Landsat-TM 
y sobre Landsat 8-OLI: 


TM: 
TM1 
TM 2 
brillo 0,3037 0,2793 0,4743 0,5585 0, 5082 0,1863 M3 
verdor |=|-0,2848 -0,2435 -0,5436 0,7243 0,0840 -0,1800 o 
humedad 0,1509 0,1973 0,3279 0,3406 -0,7112 -—0,4572 
TM5 
TM7 
OLI: 
"0,3029 0,2786 0,4733 0,5599 0,5080 0,1872] 
0 -0,2941 -0,2430 -0,5424 0,7276 0,0713 -01608 | | 9? 
Verdor OLI3 
Humedad || 0,1511 0,1973 0,3283 0,3407 0,7117 -0,4559 | | ozra 
TCT4 | [|-0,8239 0,0849 0,4393 -0,0580 0,2013 -0,2773| | OLI 
qn -0,3294 0,0557 0,1056 0,1855 -0,4349 0s808s | | 94% 
TCT6 OLIT7 
| 0,1079 0,9023 0,4119 0,0575 -0,0259 0,0252) 





Como puede verse en ambas, los coeficientes de las 
matrices de giro buscan ponderar equilibradamente 
todas las bandas para buscar la componente brillo. 
En cambio, enfatizan el contraste que existe entre el 
infrarrojo próximo y el visible, particularmente el 
rojo, para deducir la componente verdor. Finalmente, 
la humedad queda expresada contraponiendo los ca- 
nales del infrarrojo medio contra los demás. 
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La primera es obvia. El verdor se basa en la diferencia 
de reflectancia que existe entre el infrarrojo próximo 
y el rojo si la vegetación es vigorosa y el contenido 
en clorofila fotosintéticamente activa es alto. Por úl- 
timo, son muy características de la molécula de agua 
sus bandas de absorción en diversos intervalos del 
espectro en el infrarrojo medio. 


Reflectancia (%) 


407 


20- 





Comparación entre Na; 


(A/B), =2.,429 
(A/O)=2,125 ) did 


(B/C), =0,875 
(B/C), = 1,103 


] B=C 
Cocliente: 

(K/D, =0,823 

(k/h, =0,656 BC; A=C 
(K/D¿ = 0,828 
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Comparación directa de ND 
y de sus cocientes. 


Índices de vegetación 


Índices de vegetación 





En la superficie sólida del planeta solamente caben 
encontrar unos cuantos tipos de superficies genéri- 
cas: o bien se trata de suelo desnudo, o bien de vege- 
tación, o de la mezcla de ambos en cualquier propor- 
ción, o bien es agua en estado sólido o en masas lí- 
quidas. Los índices de vegetación son elementos ade- 
cuados para caracterizar el tipo de cubierta y sinteti- 
zar la información en una sola banda. Además, son 
útiles para deducir otras variables de gran trascen- 
dencia en los estudios físicos de la superficie, como 
es el caso de la emisividad. 


Cocientes de bandas 


Los cocientes son transformaciones consistentes en 
dividir celda a celda todos los ND pertenecientes a 
una banda por los de otra. El caso más general puede 
expresarse: 


(z, +4) 


zZ3=C,+C, 0 
1 


(21) 


aunque no necesariamente hay que añadir constan- 
tes aditivas. 


En ciertas ocasiones, los cocientes proporcionan ven- 
tajas en la interpretación de las imágenes. Cuando es 
difícil discriminar entre dos cubiertas cuyos ND en 
varias bandas son muy similares, el cociente entre 
ellas puede ofrecer un mejor criterio de separabili- 
dad. 


Para la misma irradiancia extraterrestre, el relieve 
origina gradientes de irradiancia local y, en conse- 
cuencia, niveles variables de radiancia que son inde- 
pendientes de la ocupación de suelo y que, además, 
afectan a todas las bandas de la imagen. Por el con- 
trario, el cociente entre bandas independiza de la irra- 
diancia la interpretación de los resultados. En conse- 
cuencia, es una transformación muy indicada para re- 
ducir el efecto radiométrico del relieve en la imagen 
final. 


Por ejemplo, en igualdad de condiciones de ilumina- 
ción, las coníferas siempre son más reflexivas que las 
caducifolias, tanto en el visible como en el infrarrojo 
próximo. En cambio, una conífera en sombra brilla 
menos que una caducifolia al sol. El cociente IRp/R 
soluciona esta cuestión. Este cociente es un índice de 
vegetación. 


Índices de vegetación 


Los índices de vegetación son transformaciones sen- 
cillas y rápidas cuyo objetivo es realzar las zonas cu- 
biertas por material fotosintéticamente activo, maxi- 
mizando la contribución espectral inherente a la cu- 
bierta vegetal, en tanto que minimizan la influencia 
de otros factores que pueden distorsionar la interpre- 
tación, como pueden ser el tipo de suelo, la irradian- 
cia solar, etc. 


Existen dos grandes grupos de índices de vegetación, 
los que se basan en transformaciones ortogonales y 
los basados en cocientes entre bandas. 
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Los índices basados en trasformaciones recurren a 
giros en el espacio espectral IRp/R. Para construir el 
Índice de Vegetación Perpendicular (PVI, Perpendi- 
cular Vegetation Index), uno de los nuevos ejes 
queda orientado según la dirección de la línea de sue- 
los, y el otro perpendicularmente a ella. El PVI mide 
la distancia espectral entre el vector de característi- 
cas de la celda estudiada y la línea de suelos: 





PVI = ¿Te (IRp, o(Rp.) | +[o(R,)-P(R,)] 
(22) 


donde p es la reflectancia en los canales IRp y R y 
los subíndices v y s se refieren a las celdas con vege- 
tación y a las de suelo desnudo. En general, las celdas 
ocupadas por suelo desnudo tendrán PVI=0, las ocu- 
padas por agua podrán tener PVI<O y las cubiertas 
por vegetación tendrán valores de PVI positivos y 
tanto mayores cuando más vigorosa sea la vegetación 
de la celda. 


Es un rasgo muy característico de la vegetación a me- 
dida que cubre el suelo una disminución progresiva 
de la reflectancia en el visible, particularmente en el 
rojo espectral debido a la absorción de la luz por 
parte de los pigmentos fotosintéticos, especialmente 
la clorofila, y un aumento simultáneo de la reflectan- 
cia en el infrarrojo próximo. Además, las líneas de 
isocobertura, para distintos tipos de atmósfera se van 
inclinando a medida que el cultivo crece y madura. 


Un índice que proporciona información acerca de la 
cantidad de masa foliar en un determinado punto es 
el índice de área foliar (LAI, Leaf Area Index) que es 
la proporción de suelo cubierto eficazmente por ho- 
jas u otro material vegetal. Conforme se desarrolla 
una cubierta vegetal, su índice foliar se incrementa 
hasta llegar a un máximo (próximo a la floración). A 
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partir de entonces decrece hasta alcanzar el punto de 
marchitez permanente. Durante el crecimiento del 
cultivo la contribución reflectora del suelo va dismi- 
nuyendo. El LAT puede superar ampliamente la uni- 
dad. 


Entre los índices basados en cocientes, los más utili- 
zados consisten en los que combinan bandas del in- 
frarrojo próximo con las del visible y se basan, como 
se ha dicho en la alta reflectancia de la vegetación en 
el IRp y la baja en el VIS, especialmente en el rojo. 
El más simple de los índices VIS-IRp es el Indice de 
Vegetación de la Razón RVI (Ratio Vegetation In- 
dex): 


Rv PE 
Pr 


(23) 


elaborado con los cocientes celda a celda de las ban- 
das IRp y R de la imagen. Si bien este índice muestra 
claramente la relación entre las relectancias de las 
bandas citadas tiene el inconveniente de no estar aco- 
tado y tender a valores muy altos cuando el denomi- 
nador es pequeño, razón por la cual no son fáciles 
las comparaciones. Por esta razón se prefiere norma- 
lizar el índice y elaborar el denominado Índice de Ve- 
getación de Diferencia Normalizada NDVI (Vormali- 
zed Difference Vegetation Index): 


Pirp Pr 
Pirp FPr 


NDVI = (24) 


La ventaja de este índice es que está acotado entre -1 
y 1. Superficies sin vegetación proporcionan valores 
de NDVI próximos a cero, tanto positivos como nega- 
tivos. Valores altos de NDVI indican celdas con alto 
contenido en vegetación vigorosa y fotosintética- 
mente activa. Valores francamente negativos de 
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Problemas en la interpretación de la cober- 
tura de una imagen si coexisten laderas bien 
y mal iluminadas. 
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Evolución espectral de una celda 
con la humedad. 


'(NDVID 


RVI 


NDVI pueden indicar superficies de agua o cubiertas 
nubosas. 


Es importante hacer notar que tanto éste como el 
resto de los índices están definidos y se calculan con 
las reflectancias, no con los niveles digitales directa- 
mente. 


Directamente relacionado con el anterior está el ín- 
dice de vegetación transformado TVI (Transformed 
Vegetation Index): 


TVI =,[[NDVI +0,5 (25) 





que evita los valores negativos del NDVI. La introduc- 
ción de la raíz cuadrada persigue asemejar la distri- 
bución de valores a una distribución normal. No hay 
grandes diferencias entre los resultados obtenidos 


con con NDVI y TVI. Se suele tomar el umbral /2 / 2 


para delimitar las áreas de vegetación. Valores supe- 
riores a esta cifra indican cobertura vegetal. 


Corrección de la influencia del suelo 


En numerosas ocasiones la vegetación no tapiza com- 
pletamente el suelo y el índice que se calcule puede 
no reflejar fielmente el estado vegetativo de la cu- 
bierta simplemente por el hecho de estar la radiancia 
contaminada por el efecto del suelo subyacente. Esto 
sucede tanto en cultivos herbáceos en hilera como en 
cultivos leñosos, pero también es frecuente en bos- 
ques del tipo mediterráneo donde la cabida cubierta 
es escasa y entre árboles y arbustos se deja ver el 
suelo desnudo. En estos casos ha de recurrirse a pro- 
cedimientos que minimicen ese efecto del suelo par- 
cialmente cubierto. 


Índices de vegetación 


La evolución espectral de una celda que paulatina- 
mente aumenta su humedad en un diagrama IRp/R 
es contradictoria en función de qué tipo de índice se 
utilice: si se utiliza el PVI, el desplazamiento se reali- 
zará según líneas paralelas a la línea de suelo, puesto 
que el contenido de material vegetal es invariable. 
Sin embargo, si se utiliza el criterio del RVI o del 
NDVI, la celda tiende hacia el origen de coordenadas. 


El índice de vegetación ajustado al suelo, SAVI (Soil 
Adjusted Vegetation Index) proporciona una solu- 
ción de compromiso entre ambas situaciones y reúne 
características comunes a los anteriores. El punto de 
convergencia de la bisectriz de cada par de isolíneas 
con la línea de suelo está en (-4, -£), y de valores 
experimentales se deduce que habitualmente (4 = £ 
= (0,25 unidades de reflectancia. Para eliminar los va- 
lores negativos en el diagrama IRp/R se traslada el 
origen a la posición (-4, - 4), obeniéndose un NDVI 
desplazado: 


Ml (Pazo +£,)-(p Est) 


NOW ensniacado = 
di Us +1,)+(Pr +£,) 

Si se particulariza este índice para los valores de in- 
tersección de la prolongación de la línea de suelos y 
la bisectriz de los pares de isolíneas (4 = £ = 0,25) y 
se multiplica por una constante que homologue su 
intervalo de existencia al del NDVI se tiene final- 
mente: 


Pirp Pr 
Pirp FPR EL 


SAVI = -(1+L) (26) 


Donde L adquiere el valor 1 para baja proporción de 
vegetación, 0,5 para una cobertura intermedia y 0,25 
para alto contenido de vegetación. 
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Otros índices de vegetación 


Existen otros muchos índices de vegetación relacio- 
nados o derivados de los anteriores, apropiados cada 
uno de ellos para circunstancias defendidas por sus 
autores respectivos. A continuación se muestra una 
mínima parte de la que los investigadores han ido 
desarrollando a lo largo de los últimos años. 


e SAVI modificado: 





2P 12p +12, +1] Ñ 81 Pr, pr] (27) 
2 


MSAVI = 
e SAVI transformado: 


al Por, =4Px =b] 
p¿+a Pig 7 ab+X (1+a*) 





TSAVI = (28) 


donde a es la pendiente de la línea de suelos, b es la 
ordenada en el origen de dicha línea y X=0,08 es un 
parámetro a determinar según el tipo de cultivo. 


e GEMI (Global Environment Monitoring): 


n(1-0,25n)-(pz -0,125) 
l-p 


GEMI = 





Donde: 


y 7 PR ) AE yn SP 12p ña 0, SPz 
Pirp sE Pr + 0,3 
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e GESAVI (Generalized SAV]I): 


Pirp 7 APR =b 
Paté 


GESAVI = (30) 


donde a y b son respectivamente la pendiente y la 
ordenada en el origen de la línea de suelos en un 
diagrama IRp/R: 


(RP ¿ero = AR ¿101 4D) 


suelo 


Z hace referencia al valor absoluto de la abscisa del 
punto de corte de las isolíneas de vegetación con la 
línea de suelo. Un valor empírico adecuado es 
Z=0,35. 


e OSAVI (Optimized SAVI): 


Pirp Pr 
Pirp FPr + E 


OSAVI = (1+Y) (31) 


donde Y es un parámetro de ajuste de valor 0,16. 
Para bajos valores de LAI, el OSAVI independiza el 
resultado del suelo subyacente muy eficazmente. Se 
utiliza sobre todo para el cálculo de la biomasa sobre 
el suelo, nivel de nitrógeno en hoja y contenido de 
clorofila, Este índice equivale al TSAVI con a=1 y 
b=0. 


e ARVI (Atmosferically Resistant Vegetation In- 
dex): 


ARVI = Pirp —[2p, Pa] 


—— ARK (32) 
Pirp + [2p, > Pa] 


en cuya expresión, pa es la reflectancia en el azul es- 
pectral. 


e CARI (Chlorophyll Absorption Ratio Index): 


Po [a -670+ Pg +bl Pro 





(33) 
Ja? +1 Pozo 


Pz00 —Psoo . 
con a b =Pss —(a-550) 


Este índice minimiza el efecto de los materiales no 
fotosintéticos en la cubierta vegetal. 


e MCARI (CARI modificado): 


2,5 (Piso = Pozo ) LS (De = Puig.) 
(28500 E 1) (Poo Pero) 0,5 


MCARI = 1,5 








(34) 


El MCART es más sensible a la concentración de clo- 
rofila en las hojas. 


e TCARI (CARI transformado): 





TCARI =3- ton Pozo ) -0,2 Ao al Pes) Po | 


670 


(35) 


Mejora la sensibilidad ante bajos contenidos en clo- 
rofila. 


e Índice de suma verde SGI (Sum Green Índex): 


Es un índice hiperespectral muy simple sensible a pe- 
queños cambios en la densidad del dosel vegetal y, 
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por tanto, válido para hacer el seguimiento de las cu- 
biertas forestales. Consiste en calcular el promedio 
de la reflectancia en el intervalo espectral compren- 
dido entre 500 y 600 nm. 


600 
Pr 
SGT=2=%"-—. (36) 
N 


Con valores 0<SGI< - 50% en porcentaje de re- 


flectancia. Para la vegetación verde el rango más co- 
mún es del 10 al 25 %. 


Índices de humedad 


Al igual que mediante cocientes se puede determinar 
el contenido en vegetación fotosintéticamente activa 
en una celda de imagen, también existen índices simi- 
lares para medir el contenido de humedad. Dos índi- 
ces representativos de ello son los siguientes, para 
cuyo cálculo se necesita conocer la reflectancia en 
ciertas bandas muy estrechas del espectro, lo que 
obliga a utilizar información hiperespectral. 


e WBI (Water Band Index): 


wB] = Po00 (37) 
Paz 


La absorción alrededor de 970 nm aumenta en rela- 
ción a la de 900 nm cuando aumenta el contenido de 
agua del dosel vegetal. Este índice se aplica al análisis 
de estrés hídrico de la vegetación, a la predicción de 
la producción, al análisis de riesgo de incendio, y en 
general en la gestión agrícola y en estudios medioam- 
bientales. Rangos habituales para la vegetación verde 
son de 0,8 a 1,2. 


Teledetección: principios y aplicaciones 


MAPA Álgebra de imágenes 


e NDWI (Vormalized difference Water Index): 


NDw] = P57—Prar (38) 
Pes7 PE Pi41 


El NDWI es sensible a los cambios en la humedad de 
la vegetación. La dispersión de la luz por el dosel ve- 
getal favorece la absorción del agua en 1241 nm. Se 
aplica al análisis de estrés del dosel forestal, al estu- 
dio del LAI áreas de vegetación densa, a los modelos 
de productividad vegetal y al riesgo de incendios. El 
NDWI está acotado entre -1 y 1, y sus valores usuales 
para la vegetación verde son de 0,1 a 0,4. 


e MSI (Moisture Stress Index): 


Ms] = Piso (39) 
Pero 


El MSI es sensible al incremento de la humedad fo- 
liar. La absorción en 1599 nm aumenta con el conte- 
nido de agua y en cambio la registrada en 819 nm 
permanece prácticamente invariable. Este índice se 
aplica al análisis de estrés vegetal, a la predicción de 
cosecha, al estudio de los riesgos de incendio y en 
estudios ambientales. La interpretación del MSI es in- 
versa en relación con otros índices: valores más altos 
indican mayor estrés hídrico. El rango puede oscilar 
entre 0 y más de 3, siendo usuales para la vegetación 
verde valores entre 0,4 y 2. 


e NDII (Vormalized Difference Infrared Index): 


NDII = Ps19 —Piós9 (37) 
Ps19 F Piso 


El NDII es sensible a los cambios en el contenido de 
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agua de la cubierta vegetal. Sus valores aumentan 
con el contenido de agua. Sus aplicaciones principa- 
les son la gestión agrícola y forestal y la detección de 
estrés hídrico de la vegetación. Al igual que el NDWI, 
sus valores están acotados entre -1 y 1, y los más co- 
munes para la vegetación verde están comprendidos 
entre 0,02 y 0,6. 


Índice de área quemada 


Se han desarrollado también índices basados en co- 
cientes para calcular la severidad de los incendios fo- 
restales. Uno representativo es el siguiente. 


e MSI (Vormalized Burn Ratio): 


NBR = Pg6s — Piosos (40) 
Psos F Progos 


Para una zona determinada, se calcula el NBR antes 
y después de un incendio forestal y se determina el 
incremento en el NBR: 


DNBR = NBR,,, — NBR 


pos 


en donde los subíndices pre y pos indican las situa- 
ciones previas y posteriores al incendio. Valores infe- 
riores a -0,25 indican una revegetación intensa pos- 
terior al incendio. Valores entre -0,25 y -0,1 indican 
baja revegetación. Valores entre -0,1 y 0,1 indican 
áreas no quemadas. Valores de 0,1 a 0,27 indican 
baja severidad del incendio. Valores de 0,27 a 0,44 
indican severidad de baja a moderada de los efectos 
del fuego. Valores de 0,44 a 0,66 indican severidad 
de moderad a alta del fuego. Finalmente, valores su- 
periores a 0,66 indican una alta severidad de los efec- 
tos del incendio. 


Si no se dispone de información hiperespectral, el ín- 
dice puede acomodarse a la anchura de banda del 
sensor TM o ETM+ de Landsat 5 o 7, o sus bandas 
equivalentes de OLI en el satélite Landsat 8: 


_TM4-TM6 


NBR = 
TM4+TM6 


Índices de vegetación MIE 


Teledetección: principios y aplicaciones 


MEA Álgebra de imágenes 


Carlos Pinilla 


